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ABSTRACT 

We show that a significant fraction of the energy of a 7-ray burst (GRB) is probably 
converted to a burst of 10 17 -10 19 eV neutrinos and multiple GeV 7's that follow the 
main GRB by ~ 10 s. If GRB's accelerate protons to ~ 10 20 eV, a suggestion that 
recently gained support from observations of GRB afterglow, then both the neutrinos 
and the gammas may be detectable. 



1. Introduction 

The widely accepted interpretation of the phenomenology of 7-ray bursts (GRB's) is that the 
observable effects are due to the dissipation of the kinetic energy of a relativistically expanding 
fireball whose primal cause is not yet known [see Meszaros (1995) and Piran (1996) for reviews]. 
The physical conditions in the dissipation region imply that protons can be Fermi accelerated to 
energies > 10 20 eV (Waxman 1995a; Vietri 1995; see Waxman 1999 for recent review). 

Adopting the conventional fireball picture, we showed previously that the prediction of an 
accompanying burst of ~ 10 14 eV neutrinos is a natural consequence (Waxman & Bahcall 1997). 
The neutrinos are produced by tt + created in interactions between fireball 7-rays and accelerated 
protons. The key relation is between the observed photon energy, e 7 , and the accelerated proton's 
energy, e p , at the photo-meson threshold of the A-resonance. In the observer's frame, 



e 7 e p = 0.2 GeV 2 T 2 , (1) 

where phenomenologically the Lorentz factors of the expanding fireball are T > 10 2 . Inserting 
a typical observed 7-ray energy of 1 MeV, we see that characteristic proton energies ~ 2 x 10 6 
GeV are required to produce neutrinos from pion decay. Typically, the neutrinos receive ~ 5% of 
the proton energy, leading to neutrinos of ~ 10 14 eV as statedQ In the standard picture, these 
neutrinos result from internal shocks within the fireball. 

In the last two years, afterglows of GRB's have been discovered in X-ray, optical, and radio 
(Costa et al. 1997; van Paradijs et al. 1997; Frail et al. 1997). These observations confirm (Waxman 



x The claim that the spectrum of neutrinos produced in interaction with burst (rather than afterglow) photons 
extends to 10 19 eV (Vietri 199S ) is due to calculational errors. 
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1997a; Wijers, Rees, & Meszaros 1997) standard model predictions (Paczyhski & Rhoads 1993; 
Katz 1994; Meszaros & Rees 1997; Vietri 1997a) of afterglows that result from the collision of the 
expanding fireball with the surrounding medium. Inserting in Eq. (|l]) a typical afterglow photon 
energy ~ 10 2 eV, we see that characteristic neutrino energies of order 10 9 GeV may be expected^. 
7-rays of similar energies are produced by 7r° decay, but because the fireball is optically thick at 
these energies the 7's probably leak out only at much lower energies, ~ 10 GeV. 

The ultra-high energy neutrinos, ~ 10 18 eV, are produced in the initial stage of the interaction 
of the fireball with its surrounding gas, which occurs over a time, T, comparable to the duration of 
the GRB itself. Optical-UV photons are radiated by electrons accelerated in shocks propagating 
backward into the ejecta. Protons are accelerated to high energy in these "reverse" shocks. The 
combination of low energy photons and high energy protons produces ultra-high energy neutrinos 
via photo-meson interactions, as indicated by Eq. (||). 

Afterglows have been detected in several cases; reverse shock emission has only been identified 
for GRB 990123 (Akerlof 1999). Both the detections and the non-detections are consistent with 
shocks occurring with typical model parameters (Sari & Piran 1999; Meszaros & Rees 1999), 
suggesting that reverse shock emission may be common. The predicted neutrino and gamma-ray 
emission depends upon parameters of the surrounding medium that can only be estimated once 
more observations of the prompt optical afterglow emission are available. 

We discuss in §2 likely plasma conditions in the collisions between the fireball and the 
surrounding medium, and in §3 the physics of how the ultra-high energy neutrinos are produced. 
We first discuss in §3.1 ultra-high energy cosmic ray (UHECR) production in GRB's in the light 
of recent afterglow observations. We show that afterglow observations provide further support 
for the model of UHECR production in GRB's, and address some criticism of the model recently 
made in the literature ( pallant fc Achterberg 1999 ). Neutrino production is then discussed in 



§3.2. The expected neutrino flux and spectrum is derived in §4 and the implications of our results 
to future experiments are discussed in §5. 



2. Plasma conditions at the reverse shocks 

We concentrate in this section on the epoch between the time the expanding fireball first 
strikes the surrounding medium and the time when the reverse shocks have erased the memory of 
the initial conditions. After this period, the expansion approaches the Blandford & McKee (1976) 
self-similar solutions. The purpose of this section is to derive the plasma parameters and the UV 
luminosity and spectrum expected for typical GRB parameters. We use the results derived here to 



2 Afterglow photons are produced over a wide range of energies, from radio to X-rays, leading to a broad neutrino 
spectrum. As we show in §4, however, the flux is dominated by neutrinos in the energy range 10 17 -10 19 eV [cf. Eq. 
which are produced in proton interactions with 10 eV-1 keV photons. 



-3 - 



estimate in the following sections the energy to which protons can be accelerated and to calculate 
the production of ultra-high energy neutrinos and multiple GeV photons. 

During self-similar expansion, the Lorentz factor of plasma at the shock front is 
Ybm = (17-E/167rnm p c 2 ) 1 / 2 r -3//2 , where E is the fireball energy and n is the surrounding gas 
number density. The characteristic time at which radiation emitted by shocked plasma at radius r 
is observed by a distant observer is t « r/4T 2 BM c (Waxman 1997b). 

The transition to self-similar expansion occurs on a time scale T (measured in the observer 
frame) comparable to the longer of the two time scales set by the initial conditions: the (observer) 
GRB duration £grb and the (observer) time Tp at which the self-similar Lorentz factor equals the 
original ejecta Lorentz factor Tj, Trm^ = Ir) = IY Since t = r/AT BM c, 

1/3 / r>. \ -8/3 



max 



(2) 



During the transition, plasma shocked by the reverse shocks expands with Lorentz factor close to 
that given by the self-similar solution, T ~ Trm^ = T), i.e. 

r~24sf— V /8 rr 3/8 , (3) 

V n ) 

while the unshocked fireball ejecta propagate with the original expansion Lorentz factor, Tj > T. 
We write Eq. (||) in terms of dimensionless parameters that are characteristically of order unity in 
models that successfully describe observed GRB phenomena. Thus E = 10 53 -E53 erg, T = lOTi s, 
n = lnocm -3 , and typically Tj ~ 300. Lorentz factors of the reverse shocks in the frames of 
the unshocked plasma are mildly relativistic, Tr — 1 ~ Ti/F^M- If the initial Lorentz factor 
is extremely large, Tj 3> 300, the transition Lorentz factor computed from Eq. (g) and Eq. (||) 
remains unchanged, T ~ 250, while the reverse shocks become highly relativistic, Tr ~ Ti/T 2> 1. 

The observed photon radiation is produced in the fireball model by synchrotron emission of 
shock-accelerated electrons. We now summarize the characteristics of the synchrotron spectrum 
that leads to ultra-high energy neutrinos and GeV photons. 

Let £ e and £b be the fractions of the thermal energy density U (in the plasma rest frame) 
that are carried, respectively, by electrons and magnetic fields. The characteristic electron 
Lorentz factor (in the plasma rest frame) is 7 m ~ S,e(^R — ^)m p /m e ~ £ e (Ti/T)mp/m e , where 
the thermal energy per proton in the shocked ejecta is ~ (Tr — l)m p c 2 . The energy density U is 
E ~ 47rr 2 'cTT 2 U , and the number of radiating electrons is iV e ~ E /TirripC 2 . The characteristic (or 
peak) energy of synchrotron photons (in the observer frame) is 

ef m - KTjI^ = V^UetA' 2 (A) 2 eV , (4) 

and the specific luminosity, L e = dL/de^', at is 

L m « (2^)-r-0iV e w 6 x 10 60 e V2 2 ^/4 T -3 /4 ^y\^ S '\ (5) 
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where £ e = 0.1£ ej _i, and = 0.01^^2- Hereafter, we denote particle energy in the observer 
frame with the super-script "ob.", and particle energy measured at the plasma frame with no 
super-script (e.g., = Te im ). 

Since the reverse shocks are typically mildly relativistic, electrons are expected to be 
accelerated in these shocks to a power law energy distribution, dN e /d^f e cx j~ p for j e > 7 m , 
with p ~ 2 (Axford, Leer, & Skadron 1977; Bell 1978; Blandford & Ostriker 1978). The specific 
luminosity extends in this case to ) 1/2 , up to photon energy 

e 7C . Here e 7C is the characteristic synchrotron frequency of electrons for which the synchrotron 
cooling time, Qirni e c/ 'ax^eB 2 , is comparable to the ejecta (rest frame) expansion time, ~ r/Tc. At 
energy e 7 > e 7C , 

e° h c » 0.3^Z 2 2 n - ^sa 172 ^ 172 keV, (6) 
the spectrum steepens to L t oc e~. 



3. UHECR and neutrino production 

3.1. UHECR production 

Protons are expected to be accelerated to high energies in mildly relativistic shocks within 
an expanding ultra-relativistic GRB wind (Waxman 1995a, Vietri 1995). Energies as high as 
€p b ' = 10 20 e° b 2 'o eV may be achieved provided the fraction of thermal energy density carried by 
magnetic fields, £b, is large enough, and provided shocks occur at large enough radii, so that 
proton energy loss by synchrotron emission does not affect acceleration (Waxman 1995a, 1999). 
The condition that needs to be satisfied by £e, 

^ > io- 2 ( e ° b 20 r/ 25 °) 2 ^52 . (7) 

where T is the wind expansion Lorentz factor and L 7 = 10 52 L 7j 52erg s -1 its 7-ray luminosity, 
is consistent with constraints imposed by afterglow observations. Afterglow observations imply 
£ e ~ 0.1 and £b ~ 0.01 [e.g. Eq. (H)]. The observed distribution of GRB redshifts, which 



suggests that most detected GRB's occur at the redshift range of 1-3 ( Krumholtz, Thorsett fe 



Harrison 19981 , [Mao fc Mo 199Sj , [Hogg fe Fruchter 1999] ), imply that the characteristic GRB 



7-ray luminosity is L 7 ~ 10 52 ergs _1 : For characteristic GRB 7-ray flux, F~ ~ 10 _6 erg/cm s in 
the BATSE 20keV-2MeV range, and adopting the cosmological parameters = 0.2, A = and 
Hq = 75km/sMpc, the luminosity for a z = 1.5 burst is L 7 ~ 10 52 ergs _1 . This result is consistent 
with the more detailed analysis of Mao & Mo (1998), who obtain a median GRB luminosity in the 
50keV-300keV range (which accounts for ~ 1/3 of the BATSE range luminosity) of ~ 10 51 ergs _1 
for Q = 1, A = and H = lOOkm/sMpc. 



The condition that needs to be satisfied to avoid proton synchrotron energy loss ( Waxmarj 
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1995a), 

r > r syn = 10 12 (r/250)- 2 (e° b 20 ) 3 cm, (8) 

is clearly satisfied in the present context, as reverse shocks arise at r 4T 2 cT « 10 17 cm > 10 12 cm. 
Thus, synchrotron losses of protons accelerated to high energy at the radii where reverse shocks 
are expected to arise are negligible. 

We note that it has recently been claimed (Gallant & Achterberg 1999) that acceleration of 
protons to ~ 10 20 eV in the highly-relativistic external shock driven by the fireball into its surround 
medium is impossible. Regardless of whether this claim is correct or not, it is not relevant to the 
model proposed in Waxman (1995a) and discussed here, in which protons are accelerated in the 
mildly-relativistic internal (reverse) shocks. 

Finally, improved constraints from afterglow observations on the energy generation rate of 
GRB's provide further support to the GRB model of UHECR production. For an open universe, 
f2 = 0.2, A = and Hq = 75km/sMpc, the GRB rate per unit volume required to account for 
the observed BATSE rate is i? G RB ~ 10" 8 Mpc- " 3 yr , assuming a constant comoving GRB rate. 
Present data does not allow to distinguish between models in which the GRB rate is evolving 
with redshift, e.g. following star formation rate, and models in which it is not evolving, since 



in both cases most detected GRB's occur at the redshfit range of 1-3 ( Hogg fc Fruchter 1999|) . 
Thus, Rqkb provides a robust estimate of the rate at z ~ 1, while the present, z = 0, rate may 
be lower by a factor of ~ 8 if strong redshift evolution is assumed. This implies that the present 
rate of 7-ray energy generation by GRB's is in the range of 10 44 erg/Mpc 3 yr to 10 45 erg/Mpc 3 yr, 
remarkably similar to the energy generation rate required to account for the observed UHECR 
flux above 10 19 eV, ~ 10 44 erg/Mpc 3 yr ( |Waxman 1995b| , [Waxman fc Bahcall 199% . 



3.2. Neutrino production 

The photon distribution in the wind rest frame is isotropic. Denoting by re 7 (e 7 )cfe 7 the 
number density (in the wind rest frame) of photons in the energy range e 7 to e 7 + de 7 , the 
fractional energy loss rate of a proton with energy e p due to pion production is 

t- 1 (e)= --^ 
K [6p) e p dt 

= — ttC / de a w (e)£(e)e / dxx~ 2 n(x), (9) 

where j p = e p /m p c 2 , a 7T (e) is the cross section for pion production for a photon with energy e in 
the proton rest frame, £(e) is the average fraction of energy lost to the pion, and eo = 0.15 GeV is 
the threshold energy. 

The photon density is related to the observed luminosity by n{x) = L e (Tx) / (A r Kr 2 cTx). 
For proton Lorentz factor eo/2e 7C <C 7 P < eo/2e 7m , photo- meson production is dominated by 
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— 1/2 

interaction with photons in the energy range e 7m < e 7 <C e 7C , where L e oc e 7 . For this photon 
spectrum, the contribution to the first integral of Eq. (^) from photons at the A resonance is 
comparable to that of photons of higher energy, and we obtain 



Here, o" poa k ~ 5 x 10 _28 cm 2 and ^ pea k — 0.2 at the resonance e = e pea k = 0.3 GeV, and 
Ae ~ 0.2 GeV is the peak width. The time available for proton energy loss by pion production is 
comparable to the expansion time as measured in the wind rest frame, ~ r/Tc. Thus, the fraction 
of energy lost by protons to pions is 



U(ef-) « 0.05 ( ( —) Tf 1 

V6 x lO^s- 1 J V250/ 1 

ob. ob. \l/2 



ob. ob. \L Z f-i -I -i 



Eq. (|11[) is valid for protons in the energy range 



4 X 1018 (2B0) ( e 7^ev)" leV < £ p b ' < 



4xl0 21 ^ 2 (6° b m , eV )~ 1 eV. (12) 



V250/ 



Such protons interact with photons in the energy range e 7m to e 7C , where the photon spectrum 

1/2 1/2 
L e oc e 7C and the number of photons above interaction threshold is oc e p . At lower energy, 

protons interact with photons of energy e 7 > e 7C , where L e oc e _1 rather then L t oc e™ 1 / 2 . At these 

energies therefore oc e° b ' . 

Since the flow is ultra-relativistic, the results given above are independent of whether the wind 
is spherically symmetric or jet-like, provided the jet opening angle is > 1/T . For a jet-like wind, 
L m is the luminosity that would have been produced by the wind if it were spherically symmetric. 



4. Neutrino spectrum and flux 

Approximately half of the energy lost by protons goes into 7r° 's and the other half to 7r + 's. 
Neutrinos are produced by the decay of tt +, s, tt + — > + — > e + + v e + + v^. The mean 
pion energy is 20% of the proton energy. This energy is roughly evenly distributed between the 
7r + decay products. Thus, approximately half the energy lost by protons of energy e p is converted 
to neutrinos with energy ~ 0.05e p . Eq. (12) implies that the spectrum of neutrinos below 
e° b ' « 10 17 (r/250) 2 (e 7 b ' keV ) _1 eV is harder by one power of the energy then the proton spectrum, 
and by half a power of the energy at higher energy. For a power law differential spectrum of 
accelerated protons n(e p ) oc e~ 2 , as expected for Fermi acceleration and which could produce the 
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observed spectrum of ultra- high energy cosmic rays (Waxman 1995b), the differential neutrino 
spectrum is n(e v ) oc e~ a with a = 1 below the break and a = 3/2 above the break. 

The energy production rate required to produce the observed flux of ultra-high energy 
cosmic-rays, assuming that the sources are cosmologically distributed, is (Waxman 1995b) 

E 2 CR dNc R /dE C R ~ 10 44 erg Mpc^r" 1 - (13) 

If GRB's are indeed the sources of ultra- high energy cosmic rays, then Eqs. (0JH) im p!y that the 
expected neutrino intensity is 

10 (ifev) <*>Vao-V*B-', (14) 

where /| 19] = 1^°% = 2) and = 1/2 for ef- > 10 17 eV and = 1 for e° b - < 10 17 eV. The fluxes 
of all neutrinos are similar and given by Eq. (0), $„ M « « « $„. Eq. (0) is obtained 



by integrating the neutrino generation rate implied by Eq. (13) and ( |ll| ) over cosmological time, 
under the assumption that the generation rate is independent of cosmic time (Waxman & Bahcall 
1999). If the GRB energy generation rate increases with redshift in a manner similar to the 
evolution of the QSO luminosity density, which exhibits the fastest known redshift evolution, the 
expected neutrino flux would be ~ 3 times that given in Eq. (£D (Waxman & Bahcall 1999). 

The neutrino flux is expected to be strongly suppressed at energy e° b ' > 10 19 eV, since protons 
are not expected to be accelerated to energy e° b ' ^> 10 20 eV. If protons are accelerated to much 
higher energy, the flux may extend to ~ 10 21 eV. At higher energy the ejecta expansion time 
rT is shorter than the pion decay time, leading to strong suppression of flux due to adiabatic 
energy loss at ef - > 10 21 Ti(r/250) 2 eV. Adiabatic energy loss of muons will suppress the and 
v e flux at e° b - > 10 19 T!(r/250) 2 eV. 



5. Discussion 

When the expanding fireball of a GRB collides with the surrounding medium, reverse shocks 
are created that give rise to observed afterglow by synchrotron radiation. The specific luminosity 
and energy spectrum describing these processes are given in Eq. (|])-Eq. (|6|). For typical values 
of the plasma parameters and fireball Lorentz factor, these equations are consistent with the 
afterglow observations. The burst GRB 990123 was especially luminous and also for this burst 
the relations Eq. (||)-Eq. (JgJ) are consistent with the observation of a reverse shock and the other 
afterglow phenomena. 

If protons are accelerated in GRB's up to energies ~ 10 20 eV, then the expected flux of 



high energy neutrinos is given by Eq. (|l4l ) and (11). Muon energy loss suppresses the Un and v e 
flux above ~ 10 19 eV, while pion energy loss suppresses the flux only above ~ 10 21 eV. Since 
protons are not expected to be accelerated to 10 20 eV (Waxman 1995a), the energy beyond 
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which the flux is suppressed will likely be determined by the maximum energy of accelerated 
protons. Measuring the maximum neutrino energy will set a lower limit to the maximum proton 
energy. The predicted flux is sensitive to the value of the Lorentz factor of the reverse shock (see 
Eq. 11), but this value is given robustly by Eq. (||) as T ~ 250. 



Will the ultra-high energy neutrinos predicted in this letter be detectable? The sensitivities 
of high-energy neutrino detectors have not been determined for ultra-high energy neutrinos whose 
time of occurrence is known to within ~ 10 s and whose direction on the sky is known accurately. 
Special techniques may enhance the detection of GRB neutrinos (see below). 

Planned 1km 3 detectors of high energy neutrinos include ICECUBE, ANTARES, NESTOR 
(Halzen 1999) and NuBE (Roy, Crawford, & Trattner 1999). Neutrinos are detected by observing 
optical Cherenkov light emitted by neutrino-induced muons. The probability P Ufl that a neutrino 
would produce a high energy muon with the currently required long path within the detector 
is P VIX rj 3 X lO^e^/lO^eV) 1 / 2 (Gaisser, Halzen & Stanev 1995; Gandhi, Quig, Reno, & 
Sarcevic 1998). Using QliD, the expected detection rate of muon neutrinos is ~ 0.06/km 2 yr (over 
2ir sr), or ~ 3 times larger if GRB's evolve like quasars. GRB neutrinos may be detectable in 
these experiments because the knowledge of neutrino direction and arrival time may relax the 
requirement for long muon path within the detector. 

Air-showers could be used to detect ultra- high energy neutrinos. The neutrino acceptance 
of the planned Auger detector, ~ 10 4 km 3 sr (Parente & Zas 1996), seems too low. The effective 
area of proposed space detectors (Linsley 1985; Takahashi 1995) may exceed ~ 10 6 km 2 at 
e v > 2 x 10 19 eV, detecting several tens of GRB correlated events per year, provided that the 
neutrino flux extends to e v > 2 x 10 19 eV. Since, however, the GRB neutrino flux is not expected 
to extend well above e v ~ 10 19 eV, and since the acceptance of space detectors decrease rapidly 
below ~ 10 19 eV, the detection rate of space detectors would depend sensitively on their low 
energy threshold. 

As explained in Waxman & Bahcall (1997), zz/s are not expected to be produced in the 
GRB. However, the strong mixing between and v T favored by Super-Kamiokande observations 
of atmospheric neutrinos indicates that the flux of and v T should be equal at Earth. This 
conclusion would not hold if the less favored alternative of ia, to Sterile occurs. 

The decay of 7r°'s produced in photo-meson interactions would lead to the production of 
~ 10 49 eV photons. For the photon luminosity and spectrum given in Eqs. (HI}), the fireball 
optical depth for pair production is higher than unity for e 7 > 10 GeV. Thus, the ultra-high 
energy photons would be degraded and will escape the fireball as multi-GeV photons. Since in 
order for GRB's to be the sources of ultra-high energy protons similar energy should be produced 
in ~ 1 MeV photons and ~ 10 20 eV protons (Waxman 1995b), the expected multi-GeV integrated 
luminosity is ~ 10% of the ~ 1 MeV integrated luminosity, i.e. ~ 10 _6 erg/cm 2 . Such multi-GeV 
emission has been detected in several GRB's on time scale of > 10 s following the GRB, and 
may be common (Dingus 1995). This is not, however, conclusive evidence for proton acceleration 
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to ultra-high energy. For the parameters adopted in this paper, inverse-Compton scattering of 
synchrotron photons may also produce the observed multi-GeV photons. More wide spectrum 
observations, of optical to > 10 GeV photons, are required to determine whether the observed 
multi-GeV emission on ~ 10 s is due to inverse-Compton or 7r° decay. 

We note here that Multi-GeV photon production by synchrotron emission of ~ 10 20 eV 
protons accelerated at the highly-relativistic external shock driven by the fireball into its surround 
medium has been discussed in Vietri (1997b) and Bottcehr & Dermer (1998). However, protons 
are not likely to be accelerated to such energy at the external shock (Gallant & Achterberg 1999), 
and, moreover, even if acceleration is possible, the fraction of proton energy lost by synchrotron 
emission at the radii where external shocks occur is <C 1, [see discussion following Eq. (Js|) ] and 
hence the expected flux is much smaller than the inverse-Compton or tt° decay flux. 

JNB was supported in part by NSF PHY95-13835. EW was supported in part by BSF Grant 
9800343, AEC Grant 38/99 and MINERVA Grant. 
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